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ABSTRACT 

We present an analytical theory of thermonuclear X-ray burst atmosphere 
structure. Newtonian gravity and diffusion approximation are assumed. Hydro- 
dynamic and thermodynamic profiles are obtained as a numerical solution of the 
Cauchy problem for the first-order ordinary differential equation. We further 
elaborate a combined approach to the radiative transfer problem which yields 
Q\ \ the spectrum of the expansion stage of X-ray bursts in analytical form where 

Comptonization and free-free absorption-emission processes are accounted for 
>■ ; and r ~ r~ 2 opacity dependence is assumed. Relaxation method on an energy 

; opacity grid is used to simulate radiative diffusion process in order to match ana- 

lytical form of spectrum, which contains free parameter, to energy axis. Numeri- 
cal and analytical results show high similarity. All spectra consist of a power-law 
soft component and diluted black-body hard tail. We derive simple approxi- 
mation formulae usable for mass-radius determination by observational spectra 
Oh" fitting. 
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^ ' 1. Introduction 



First discovered by Grindlay et al. (1975), strong X-ray bursts are believed to occur 
due to thermonuclear explosions in the bottom helium-reach layers of the atmosphere accu- 
mulated by a neutron star during the accretion process in close binary system. Since then 
dozens of burster-type X-ray sources were found. One of the distinctive feature of Type I 
X-ray bursts is the sudden and abrupt (~ 1 s) luminosity increase (expansion stage) fol- 
lowed by exponential decay (contraction stage). Energy released in X-ray radiation during 
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the first seconds greatly exceeds the Eddington limit for layers above the helium burning 
zone which are no longer dynamically stable. Super-critically irradiated shells of atmosphere 
start to move outward, producing an expanding wind-like envelope. The average lifetime of 
an X-ray bursts is sufficient for steady-state regime of mass loss to be established when the 
local luminosity throughout the most of the atmosphere is equal or slightly greater than the 
Eddington limit. 

During the last two decades the problem of determining properties of radiatively driven 
winds during X-ray bursts was subjected to extensive theoretical and numerical studies. 
Various theories were put forward with gradually increasing level of accuracy of the problem 
description, but only a few approaches addressed the case of considerably expanded photo- 
sphere under influence of near-Eddington luminosities (London et al. 1986; Ebisuzaki 1987; 
Lapidus 1991; Titarchuk 1994). See Lewin, Paradijs & Taam (1993) for a detailed review 
of X-ray burst study during 80's and the beginning of 90's. 

Similarly to the problem of accretion flows, notion of the existence of sonic point in 
continuous flow became a natural starting point in the analysis of wind flows from stellar 
objects. Ebisuzaki al. (1983), hereafter EHS, investigated the structure of the envelopes with 
steady-state mass outflow and pointed out the higher Eddington luminosity in the inner shells 
due to the prevalent higher temperatures and correspondingly lower Compton scattering 
opacities. They showed that the product of opacity and luminosity remains almost constant 
throughout the atmosphere which is the key assumption of the model. The existence of wind- 
like solutions for critically irradiated atmospheres was proved. Titarchuk (1994), hereafter 
T94, studied analytically spectral shapes of the expansion and contraction stages of bursts. 
He showed how EHS's approach to hydrodynamic problem can be greatly simplified with 
the sonic point condition properly calculated and tied with conditions at the bottom of the 
envelope. Haberl & Titarchuk (1995) applied the T94 model to extract the neutron star 
mass-radius relations from the observed burst spectra in 4U 1820-30 and 4U 1705-44. 

Nobili et al. (1994), hereafter NTL, adopted a high accuracy numerical approach to 
the problem of X-ray burster atmosphere structure based on the moment formalism (Thorne 
1981; Nobili, Turolla & Zamperi. 1991). They integrated a self-consistent system of frequency- 
independent, relativistic, hydrodynamical and radiative transfer equations over the whole 
atmosphere including the inner dense helium-burning shells. Three important characteristic 
of X-ray burst outflow were obtained in this work: the helium-burning zone temperature 
was maintained approximately at the level of 3 x 10 9 K, the temperature of the photosphere 
was shown to depart appreciably from the electron temperature and to stay constant at the 
outer shells, and the existence of the maximum and the minimum values of the mass loss 
rate was found. 
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One of the goals of these studies was to provide the algorithm of determination of 
the compact object characteristics by analyzing observational data. With the advent of 
high spectral and time resolution observational instruments (such as Chandra, RXTE, USA, 
XMM-Newton missions) the task of obtaining a suitable tool for fitting the energy spectra 
became extremely important. Despite numerous earlier studies of X-ray burst observations, 
recent developments have shown a growing interest of the astrophysical community in this 
area (Strohmayer & Brown 2001; Kuulkers et al. 2001). 

Obviously, the problem of radiative transfer in relativistically moving media is very 
complicated one and under rigorous consideration it must be solved numerically. In this 
paper we develop an alternative approach which allows both numerical and analytical so- 
lutions and successfully accounts for all crucial physical processes involved. We show how 
this problem under some appropriate approximations yields the spectrum of radiation from 
spherically symmetric outflows in an analytical form. We concentrate on the case of extended 
atmosphere with inverse cubic power dependence of the number density on radius, which is 
more appropriate for the expansion stage but can also be employed for description of the 
contraction as a sequence of models with decreasing mass-loss rate. 

We represent a numerical approach to the problem which then provides the validation of 
our analytical description. We adopt the general approach formulated in EHS and developed 
in T94. The problem of determining profiles of thermodynamic variables of steady-state 
radiatively driven outflow was solved in T94. The problem is reduced to the form of a 
first order differential equation, which allows easy and precise numerical solution. For a 
completeness we present this method in Section 2. Using atmospheric profiles obtained for 
different neutron star configurations we solve the problem of radiative transfer by relaxation 
method on an energy-opacity logarithmic grid. We perform temperature profile correction 
by applying temperature equation to the obtained spectral profiles. The basic formulae 
are given in Section 3.1. Then the analytical description of the problem is represented 
in detail. The analytic solution of radiative transfer equation on the atmospheric profile 
t ~ r~ 2 is presented in T94. Here we review the solution by carrying out the integration 
without introducing any approximations. In Section 4 we compare and match our analytical 
and numerical results to describe the behavior of free parameter. We finalize our work by 
examining the properties of our analytic solution, combine it with the results of Section 4 
and construct the final formula for fitting the spectra in Section 5. The discussion of our 
method along with some other important issues concerning the problem being solved are 
presented in Section 6. Conclusions follow in the last Section. 
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2. Hydrodynamics 

As we already mentioned, that the calculation of X-ray burst spectra can be treated as 
a steady state problem. To justify this assumption one has to compare the characteristic 
times of phenomena considered. Time scale for the photosphere to collapse can be estimated 
as follows 



tcoii= / , where v coU ^\ (1 - I). 

Jr s Vcou V r 

Here r s denotes the sonic point radius, which is adopted as the outer boundary of the pho- 
tosphere throughout this paper. Dimensionless luminosity is I = L/LEdd, where Eddington 
luminosity is given by 

AixcGM ns 

J^Edd — • UJ 

K 

Opacity k is expressed by the Compton scattering opacity with Klein-Nishina correction by 
(Paczynski 1983) 

«= *° „ (2) 

kq = 0.2(2 — Y He ) cmV 1 with Y He being the helium abundance, and a = 2.2 x lO -9 !^ -1 . 
It is exactly this temperature dependence of the opacity that is responsible for the excessive 
radiation flux, which appears to be super-Eddington to the outer less hot layers of the 
atmosphere. In the framework of strong X-ray bursts the following condition are usually 
satisfied: r s > 10 3 km >> r ph , I ~ 0.99. Putting m = M ns /M Sun ~ 1 results in a time 
of collapse of the order of several seconds, the observed time which a Type I X-ray burst 
usually lasts. For evaluation of the time for photons to diffuse through the photosphere, we 
note that a number of scattering events is iV m [see, for example, Rybicki & Lightman 
(1979)], where r p h is the total opacity of the photosphere, which is ~ 10. The time for a 
photon to escape is 

r 2 r 

t eS c ~ — ~ — Tph ~ 0.1 S. 

aTn e c c 

This indicates that the hydrodynamic structure develops at least ten times slower than the 
photons diffusing time through the photosphere. Although these time scales can become 
comparable in cases of greatly extended atmospheres, generally steady-state approximation 
is acceptable. 
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2.1. Basic equations for radiatively driven outflow 

The problem of mass loss as a result of radiatively driven wind was formulated by EHS. 
For convenience of the reader we summarize all equations important for the derivations in 
the following sections and refer the reader to EHS paper for details. The system of equations 
describing steady-state outflow in spherical symmetry consists of a well known Euler (radial 
momentum conservation) equation 

v dv + GM + ldP _ Q 
dr r 2 p dr 



the mass conservation law 



-^-Unr 2 pv) = 0, (4) 
dr 



the averaged radiation transport equation in the diffusion approximation 



and the entropy equation 



_ mracr^dT 
KU ~ 3~p (5) 



ds 1 dL r . . 

vT T- + 1 — 2 — = °' 6 
dr An\r A p dr 



where P, p, T, s and L r are, respectively, the pressure, the density, the temperature, the 
specific entropy, and the diffusive energy flux flowing through a shell at r. 

The outflowing gas is taken to be ideal. Dimensionless coordinate y, which is the ratio 
of the radiation pressure P r , to the gas pressure P g , is introduced by 

y_Pr_ pm p aT 3 /3 ^ 

Pg k P 

where p, — 4/(8 — 5Yn e ) and m p are the mean molecular weight and the mass of proton, 
respectively. Then the other thermodynamic quantities are expressed in terms of y and T as 



3k y ' 
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' [4y + lny-(3/2)lnT], (10) 



h = —(4y + 5/2)T, (11) 

where h is the specific enthalpy. 

Integrals of equation (4) and (3) give energy flow and mass flux correspondingly 

Anr 2 pv = $, (12) 

*2M. + h ) 9+Lr = 1l . ( i3) 



To make two more integrations, which can not be performed analytically, the constancy 
of K,L r , which stands for the integral of equation (5), over the relevant layers is assumed. 
In EHS this assumption is confirmed numerically. We can also justified by the following 
consideration. At the near-Eddington regime the radiation pressure aT 4 /3 is much greater 
than the pressure of gas almost everywhere except for the innermost layers adjacent to the 
helium-burning zone. Neglecting the gas pressure in equation (3) and multiplying it by — r 2 
we get 

^=GM + r 2 ^. (14) 
47rc dr 

Here we moved first two terms of (3) to the right hand side and used equation (5) to express 
third term by kL t . For inner and intermediate layers of the atmosphere the last term in 
(14) is negligible and equation reduces to kL v = KqLq. This term can become considerably 
large for the outermost layers where L r must exceed Lem- This also is in agreement with 
observations of X-ray bursts from which super-Eddington luminosities are inferred. For 
the sake of analytical consideration we consider nL r to be constant throughout the whole 
atmosphere, and the third integral is 

kL t = k L = const. (15) 

Replacing L r of equation (6) with equation (15), the fourth integral is obtained as 



const. 



(16) 
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Boundary conditions need to be imposed at the bottom and the outer boundary to 
determine the four integration constants $, L , S and to obtain a specific solution. 

At the bottom of the atmosphere close to the helium-burning zone there should be a 
point where gas and radiation pressure are equal. As another important numerical result 
EHS showed that near the neutron star surface the temperature and radius profiles level off 
with respect to y, so there is always a point where 

r = r b , T = T b , y= 1, (17) 

and r b is well approximated by the radius of the neutron star R ns . However, T b can not be 
considered as a real temperature of helium-burning shell at the bottom of the star surface 
because thermonuclear processes are not included in the model. Rigorous account of helium- 
burning NTL shows that temperature of burning shells vary in small range of values. 

To obtain the outer boundary condition the concept of sonic point is used. For the 
solution to be steady-state and to have finite terminal velocity it should pass sonic point 
where 

2 GM ns fdPA ( k 



»'=^={at)r{^) Y - T -' (18) 

_ A + 4(l + y,)(l + 4y a ) 

A + 3(l + % s ) ' li9j 

where A is a quantity related to the ratio of the energy flux to the mass flux (see formula 
22, below). In EHS this formula contains a typo. We give a proper derivation of this form 
for Y s in Appendix C. 



2.2. ODE solution of the hydrodynamic problem 

T94 has shown how the treatment of the hydrodynamic problem can be reduced to a 
Cauchy problem with the boundary condition determined at the sonic point. This treatment 
provides a high-accuracy method of obtaining the hydrodynamic solution. The crucial point 
is to relate the position of the sonic point with the values of the velocity and the thermody- 
namic quantities before solving the set of appropriate hydrodynamic equations. The profile 
of the expanded envelope is then obtained as a result of the integration of a single first-order 
ordinary differential equation (ODE) from the sonic point inward up to the neutron star 
surface. For completeness we present the details of this approach. 
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At the bottom of the atmosphere the potential energy per unit mass of the gas, GM/r, 
is significantly greater than the kinetic energy, v 2 /2, and enthalpy. Therefore, by ignoring 
these terms in equation (13), we obtain the value of the mass flux 



Rn 



GM„ 



-aT h L% 



b-^Edd- 



(20) 



The inner boundary condition (17), the integral (16), and equation (10) for entropy can be 
used to find the temperature distribution with respect to y 



T = T b y- 1 / x exp 



4(y-i) 

A 



and 



apirrip 3 
A = ^$ °~2- 

The condition at the sonic point (18) allows us to find the constant \1> 

T , -r r , s 3 GM 

= h s <$> + L r {r s ) - $. 

4 r. s 



(21) 



(22) 



(23) 



Combining the mass and energy conservation laws (12) and (13), the specific enthalpy and 
density equations (9) and (11) and eliminating the radial coordinate r between equations 
(12) and (13) yield the following dependence of the velocity derivative v' with respect to y 



1 + 3 



1 + % 
A 



1 _ 75 2 rT(8-5Y He )(l + 4y)(l + aT) 
y • \r bfi T b ^/^-v 2 /2-h + GM ns /r) 



(24) 



Derivation of equation (24) is given in Appendix D. 



By imposing boundary conditions at the bottom of the extended envelope (at the neu- 
tron star surface) and at the sonic point, we can determine the four integration constants 
necessary to obtain a specific solution. One can note the obvious fact that the bottom of the 
envelope can not serve as a starting point of integration of equation (24) as long as v b = 0, 
which introduces uncertainty. Fortunately, we can calculate parameters at the sonic point in 
the framework of our problem description by solving a nonlinear algebraic equation, which 
involves only y s , the ratio of the radiation pressure P r to the gas pressure P g at that point. 
Substitution of the radial coordinate r s and velocity v s from the definition of the sonic point 
position (18), and the sonic point density, p s from equation (9), we find 



GM ns \im v 
2kY s T s ' 



k 



-Y X T, 



Ps 



ajj,m p 



T 3 

s 



3k y s 
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and the expression for the temperature given in equation (21) into the mass conservation 
law (12), after some algebra give an equation for the value of y s 



(2-Y He )m 2 V /3 T u 



0.149 {%-bY He )s/*Y a y\' 



l/A+2/3 



+ 1. (25) 



Here and T^g are the neutron star surface radius and temperature in units of 10 6 cm and 
10 9 Kelvin respectively. Since Y s is expressed in terms of y s (eq. 19), equation (25) can be 
solved to determine the value of y s . Knowledge of y s can then, by substitution in equation 
(21), yield the value of the temperature at the sonic point T s and then v s from equation (18). 
It is now possible to relate v s to T s , T b , r s , r b , thus obtaining the analytical expression for 
the various dynamical quantities at the sonic point in terms of the values of the parameters 
associated with the boundary conditions. To obtain the solution of the hydrodynamical 
problem for a particular set of input parameters we use a standard Matlab/Octave package 
function minimizators and ODE solvers. 



3. Radiative Transfer Problem 

The radiation field of X-ray burst atmosphere may be described by the diffusion equa- 
tion, written in spherical geometry, with the Kompaneets's energy operator (see T94): 



1 ( d 2 J v 2 dJ v \ _ _ B j_ 



3 V Or 2 ro>T Or J cut 



/•••/; a ( dj v t 

x o^-[x —-3J u + —J u \, (26) 



m e c 2 dx \ dx T 

where xq = hv/kTo is a dimensionless frequency, To being the effective temperature; en// and 
oit = o"T^ e are the coefficients of free-free absorption and Thompson scattering, respectively, 
whose ratio is given by (Rybicki & Lightman 1979) 

^ = 1.23^(1 - W2) 7/8 *(*o) (I)' 72 ( 27 ) 

with 

*(x ) = ~ 9{XoT ° 3 /T) (1 - e^' T ) . 



X'n 
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Here g(x) is the Gaunt factor (Greene 1959) 

V3 



g(x) = —e*/ 2 K (x/2), 



Kq(x) is the Macdonald function, and gu denotes the free-fall acceleration onto the neutron 
star surface, in units of 10 14 cm s _1 . 



We combine equation (26) with the outer boundary condition of zero energy inflow 

'dj v _ 3 



dr 



= 



r=0 



(28) 



and the condition of equilibrium blackbody spectrum at the bottom of the photosphere, 
which is represented in a dimensionless form as 



B„ 



Xr 



exp(x Q T Q /T) -1" 



(29) 



We will make use of the temperature equation, which is obtained by integration of (26) 
over frequency. The opacity operator vanishes as a result of the total flux conservation with 
respect to optical depth, leaving us with 



/,•:/; ( A • 



T, 



J u dx - — I x J u dx 
o 1 Jo 



= 1.23 pgl{\l- W2) 7/8 



T 



1/2 



J u ty(x )dx 



2\/3 T 

7T Tn 



(30) 



In the condition of the extended photosphere of X-ray bursts, the density usually is very low 
and the left-hand side of the last equation can be neglected reducing the last equation to 
the formula for temperature 



T 
To 



oo poo 

x J u dx 1 J J u dx } 



(31) 



We will use the last relation to produce a corrected temperature profile for the photosphere 
where it departs significantly from that given by the hydrodynamic solution. 



3.1. Analytic Description of Radiative Diffusion 



Hydrodynamic profiles calculated in section 2 show that during the expansion stage in 
the vicinity of the sonic point v ~ v s (r/r s ). Considering this relation to be true throughout 
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the entire envelope, according to the mass conservation law, we can write 

Y He \ $ GM$ 
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inrripvr 2 87iVgm p 



1 - ^ ) r- 3 



(32) 



where <£> is the mass loss rate and v s is the velocity of gas at the sonic point. In this case 
opacity can be expressed as 

Cut 



T = C 



dr 



2r 2 



Noting that in this case 



t = 



IT' 



we can rewrite the radiation transfer equation in the form 



d 1 d J v 3 a// 



J v — B v ) 



3kT P 



dr r dr t a T 
The boundary conditions are given by 

Mr=r th = B u {r th ) 
at the inner boundary of photosphere, and 



m e c 2 r 



L V (J U ). 



3 jo 



or 



=o AnR 2 



(33) 



(34) 



(35) 



(36) 



(37) 



at the sonic surface. The ratio Oiff/ax can be written in the form 



a ff 



= 1.23 1 



Y, 



He 



-5/8 



2m T 

7/2 



3/2 



3 \ 1/2 



xa -r/. gW(Lz£3 (I) T 3 /2 = DS(l)T 3/ 2> 

where x = hv/kT e) and ^(x) = g(a;)(l — e 



(38) 



/ x 3 . 



Stated in this way the problem of radiative transfer allows an analytical approach. The 
solution of the radiative transfer equation (35) is 

t 8/7 



J(t,x) = B v - i 



' 2 4/7 r (nj 

where K p is the modified Bessel function of the first kind, and 



Hh) 



(39) 



t=-y/3DV(x) r 7/4 . 
Details of derivation of this formula are given in Appendix A. 



(40) 
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3.2. Evaluation of r th and T c 

The next step is to find the color temperature and to determine the thermalization depth 
Tth where the boundary condition (36) is valid. For saturated Comptonization, the occupa- 
tion number behaves in accordance with Bose-Einstein photon distribution n = (e M+x — 1) _1 
which might be described as a diluted blackbody spectrum or as a diluted Wien distribution. 

At first we evaluate the color temperature assuming a blackbody spectral shape. We 
look for the solution of the form 

n(r,x) = ^\- (41) 
e x — 1 

The solution, which is described in detail in Appendix B, gives for R(r) 

7/4 N 



2 3 / 7 T 

flM = 1 " f(I) ^ 



T 



Tth 



(42) 



As long as R(t) = 1 for r > T t h, there is radiation equilibrium for optical depths deeper 
than the photospheric envelope. The temperature equation in the zone < r < r t h reads 



T \ 4 2H 



T J i? 2 



where H = (47ri? 2 ts 7r 5 /15)/167r 2 and r R (< 1) is the optical depth coordinate at the outer 
boundary of the expanded atmosphere, r = R (see eq. 35). This equation can be rewritten 
as follows 

Tj 2T n ,R( T ) [ ' 

Neglecting r R with respect to r and making use of Taylor expansion (B6) of R(t) we get a 
constant value of the temperature 

T \ 4 3 2 8 / 7 r (^) t 2 t 2 

t) ^ H-RW 1 f^ - 3rM ^ (44) 

Using the notation 



formula (38) becomes 



p = 2g(x*) w In ( — ) (45) 



frp \ 7/2 

D = dA4) (46) 



T J 
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where 



D n = 1.23 1 



Y He \~^ [2m p \^ ( V /2 _ 7/8 



while we can write for r t h 





"4 /"-"" 


-4/7 


6 2^ 


-2/7 y 


r 6 2 " 


Tth = 


.7^. 




49 F 




49 P 



-2/7 



-2/7 



(47) 



Substituting it into (44) we get for T/T 



T 

— = 0.596 
T 



49^ 



2/7 --2/7 



1/2 



(48) 



Assuming the same electron number density as in (32) we than express opacity at the neutron 



star surface r ns in the form 



T ns = I 1 



2/W VStt^J ' 



(49) 



If we use the dependence of input parameters g u and $ on m, r fei6 ,T fe 9 and Y He , the next 
useful equations for the color ratio T/Tq, color temperature kT, and thermalization depth 
T t h are found: 



T _ 0.191(2 -Y He )Wr%vl¥ 4 

T " m 3/28 T 6 ^ 14 p 4 / 7 

fcT = 0.4mVV- 6 5/14 T- 5 /Xf 4 (2-^e)- 3/1 V 4/7 keV , 
r to = 90.5m 2 /V- 6 3 / 14 T b - 3 / 14 </ 8 14 (2 - Y He f^'\ 



(50) 

(51) 
(52) 



Here v s $ is the sonic point velocity in units of 10 8 cm/s. These relations present the final 
results of our analytical approach. There is still a lack of completeness due to the presence 
of p and v s in the left-hand sides of this system of equations. Parameter p and sonic point 
velocity are not independent parameters of the problem, but at this point they can not be 
inferred from further analytical consideration. Fortunately, these quantities can be quite 
well approximated by a power dependence from m, r^T^g and Yn e , which is done in the 
next chapter. 
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Fig. 1. — Examples of spectra (left) and temperature profiles (right) obtained for model 
with m = 1.5, R ns = 13.5 km, Yue = 1- On the left solid line represents the analytical 
solution, dashed line indicates the diluted blackbody level, '+' - results of relaxation method 
simulation. On the right solid line is the temperature profile obtained from the initial hydro- 
dynamical solution, dashed line is the corrected profile (see text), dotted line is analytically 
calculated color temperature level. 




Fig. 2. — Same as Fig.l but for the model with m = 2.6, R ns = 14 km, Y He = 1. 
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4. Numerical Results and Comparison with Analytical Description of 

Radiative Transfer Problem 



To confirm the validity of our analytical approach and to examine the behavior of p 
and v s in dependence of different input parameters of the problem, we perform numerical 
modeling of the steady-state radiative transfer process. The whole procedure consisted of 
three steps. 

First, for particular model of neutron star, i. e. for a given mass and radius, we obtain a 
set of model atmospheres for a chosen set of bottom temperatures. These solutions provide 
us with runs of thermodynamical and hydrodynamical profiles, sonic point characteristics, 
masses of the extended envelopes and their loss rates. Second, we solve radiative transfer 
equation (26) on each model atmosphere obtained with the relaxation method (e.g. Press et 
al. 1992) on an energy-opacity grid using logarithmic scale on both dimensions. The energy 
range included 500 grid points. The number of grid points in opacity varied between 100 and 
300. The opacity domain included the range r s < r < r max , where r s is opacity at the sonic 
point and r max was taken large enough to meet safely inequality r max > T t h- We used the 
mixed outer boundary condition (28). Spectrum at the inner bottom of the photosphere was 
taken a pure black body B u . Numerical calculations of frequency-dependent radiation field 
consisted of two runs of our relaxation code. The first run was performed on temperature 
continuum, obtained from the hydrodynamical solution (see section 2.5). Then we calculated 
a spectral temperature profile using formulae (31) which exhibits a quite expected behavior. 
At some region this corrected profile departs from initial temperature profile and levels off 
at some constant value in absolute agreement with analytic result of section 3.2. It is also 
in a qualitative agreement with the NTL self-consisted calculation of radiation driven wind 
structure of X-ray burster. A second run is performed on the corrected profile to get more 
reliable spectrum shape. At the final third step we compared analytical and numerical 
solutions. The sonic point provided a natural position to match numerical and analytical 
solutions. Combining the sonic point parameters, calculated through the solution of equation 
(25) and using relation (34) we get for the opacity at the sonic point 



We calculated and and plotted fluxes given by both methods at the sonic point. A 
particular value of parameter p for analytic model was obtained by matching value of kT 
and corrected level of numerically achieved value of photospheric temperature. 

We obtained results for approximately 150 different sets of values T b , R ns , M ns and Y He . 
Examples of numerical calculations of spectra for different neutron star models and fitting 




(53) 
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them with analytical shapes are presented in Figures 1 and 2. Analytical and numerical 
shapes match quite well in the wide range of neutron star surface temperatures and both 
show two distinctive features of outgoing spectrum of expansion stage: diluted black-body 
like high-frequency component and power law soft excess at the lower part. Dependence of 
the sonic point opacity presented by (53) describes correctly the dilution process indicating 
the assumption of atmosphere structure adopted at the analytical model is correct. 

Tables 1 and 2, which summarize results for two different neutron star models, are given 
in order to compare our results with more rigorous calculations (NTL). Taking mass-loss rate 
as an input parameter, NTL obtained profiles of different quantities throughout the whole 
atmosphere. They argue that the temperature of burning shell is maintained around 3 x 10 9 
K for all models. The temperature of photons departs appreciably from the temperature 
of ambient matter above photospheric radius and stays practically constant indicating that 
radiation becomes essentially decoupled from expanded media. We changed the bottom 
temperature in a wide range of values and inferred the mass-loss rate, the mass of envelope, 
etc. 

Our results are in qualitative agreement with NTL. The crucial physical parameters 
which define the main spectral signatures are the photospheric radius r p h and its temperature 
kT. Runs of the atmospheric profiles obtained by both approaches are quite similar, although 
T p h in NTL results is usually 15% ~ 25% greater than in our models. This difference is 
explainable. We match isothermal levels given by numerical and analytical calculations and 
define the obtained value as a photospheric temperature. This is the lowest estimation, 
because the temperature profile starts to grow before the bottom of the photosphere is 
reached. NTL define T p h as a matter temperature at r p h- A temperature level calculated at 
the thermalization depth r th should compensate the considered difference. The difference in 
density profiles, which can achieve a factor of two, will affect the spectrum only in the soft 
part (< 0.2 keV) where the normalization of the power law component can be changed. This 
fact does not diminish the validity of our results. The soft component of the spectrum can be 
represented as an independent fitting shape with a normalization included as an additional 
fitting parameter. This matter is not crucial at the moment due to the restricted spectral 
capabilities of current X-ray observing facilities. One can also notice a quick decrease of the 
envelope mass, and point out a wide variation of T b . This discrepancy can be explained by 
differences in model formulations. Specifically, NTL included helium-burning shells into the 
model and put the inner boundary condition on the "real" neutron star surface while our 
model stops where radiation and gas pressures are equal (y = 1), which is close but still 
outside of the helium-burning shell. In our approach, part of the bottom of the atmosphere 
is left out. In fact, the lower the mass-loss rate, the greater the portion of mass missing 
beyond the point where y — 1. This is clearly seen from the tables. The temperature at the 
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bottom can be considered as an "effective" instead of the real temperature of helium-burning 
zone. 

As we have already point out that one needs to know the dependencies of v s and p on 
input parameters to complete the analytical description and thus to employ these results 
to the fitting of observational X-ray spectra. Analysis of v s and p runs show that \ogv s 
and logp are linear functions of logT^, logi? ns , logM ns and log (2 — Yue)- We combine all 
experiments and fit i> S;8 and p with a model cons x T fe a 9 r^ 6 m 7 (2 — Y He ) v by the least squares 
method to get 

p = 7.69T b -°- 84 r- 6 °- V 69 (2 - Y He )~^ (54) 

v s>8 = 5.46T^- 71 r M °- V 63 (2 - Y He )-^ (55) 

with maximum errors of parameters less than 1%. The ranges of parameters included in 
fitting are 0.3-7.0 for T^g, 0.6-2.0 for r^, 0.8-2.7 for m and 0.3-1.0 for Yn e - These results 
can be used to substitute p and v St $ in equations (50)-(52). Now we have consistent system 
of equations, which should yield X-ray spectrum of burster in the form of function of only 
input physical parameters, i. e. neutron star mass, radius, surface temperature and elemental 
abundance. 



5. Final form of the profile for spectral fitting 

The fact that spectra obtained are a black-body like almost everywhere except for small 
values of energies allows us to proceed with simplification of the formula (39). First we note 
that due to (47) and smallness of x 

4 , — — 7/4 J2V(x) J\n(2.3h/x) , N 

t th = -v^prf = 2^ U. ~ 2 ^ — 56 

7 p px 

for the soft part of spectrum. Here x = hu/kT, fy(x) and D are defined in formulae (27) 
and (46) correspondingly. Because t t h is large for small values x we can use approximation 
of the modified Bessel function of the second kind for large arguments 



and rewrite equation (39) as follows 

\ 2 

T \ 



r ( 3 / 7 ) ^8/7 + 4 ~ 1/14 C " 22 



r(ll/7) 7T(ll/7) 



-18- 



= Bv (j-^J 2 [2.32 z s ' 7 + 0.64 z 1 ' 1 ^] , (57) 



where 



v/ln(2.35/a;) 
z = . 

px 

Here we rewrite the dilution factor in terms of opacity using relation (40). Clearly, the second 
term in the parenthesis of formula (57) is significant only where z becomes small (x becomes 
large) and the spectrum shape "adjusts" to the black-body component. In turn, the first 
term of equation (57) represents the power law component of the lower part of spectrum 
with the slope 6/7, which can be shown by simple similarity (see also T94) 

X s ' 7 

Another important advantage of this term is that it vanishes for large values of x. This fact 
gives us opportunity to construct convenient and accurate formula for observational spectra 
fitting. We drop the second term in equation (57) and adjust to the diluted black-body 
shape by means of quadratic power combination as follows 



J{t,x) = B v Q-) 2 [1 + 5.34 z™' 7 ] 1/2 . 



(58) 



Comparison of shapes given by formula (58) with exact solution (39) shows that they deviate 
from each other by less than 2% which is more than acceptable in contemporary astrophysical 
observational data analysis. Using the explicit form of z and the form of outgoing flux 
equation (58) can be rewritten in the form 



47T dj u 87T n , . 

'th 



F ' = TW = T B "W< 1 + SM 



ln(2.35/x) 
p 2 x 2 



8/7 ^ V2 



(59) 



Equation (53) yields useful relationship for the dilution coefficient in the manner similar to 
(50)-(52) 



(60) 



87r r s _ 5.07 x io-vyryy^ 

Substituting results of parameter fitting (54)- (55), we get 

= 3.31 x 10- 3 r M - 36 T 6 - - 29 m-°- 17 (2 - Y^r ' 58 , (61) 

T th 
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V = 7.69T 6 / 8 V V 69 (2 " Y He y^. (54) 

Here, again, r^, T^g and m are the neutron star radius, surface temperature and mass in 
units of 10 6 cm, 10 9 Kelvin and solar mass respectively. Now we have provided our spectrum 
profile (59) with expressions for parameter p (54) and the dilution factor (61). These three 
formulae constitute the final analytical results of this paper. 



6. Discussion 

The model and derivations presented above assume that plasma consists of fully ionized 
hydrogen and helium. In reality, this assumption can be too simplistic. For instance, in the 
case of the recently discovered super-burst (Strohmayer & Brown 2001), a sufficient fraction 
of material should be represented by heavier elements. These long and powerful bursts are 
also considered to be due to the nuclear runaway burning in the carbon "ocean" under the 
neutron star surface. In this section we discuss how our model can be adjusted for study of 
this phenomenon. The approach as a whole does not change, but some formulae have to be 
modified in order to account for the different plasma composition. 

First, we note that for plasma which consists of a single ionized element, we have for 
the mean molecular weight 

A 

p = 



1 + Z' 

and for the electron number density 

n e = -f—Z, 

Am p 

where A and Z are the atomic weight and the atomic number of the corresponding element. 
In the general case of heterogenous elements, each represented by weight abundance Y^ we 
write 

" = E«(ilw (62) 

p 

n, 



In the hydrodynamic part of this study, these modifications will affect only the form of the 
terms and factors containing Yue- In the radiation transfer section, the form of aff/ar will 
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require more careful treatment. According to Rybicki & Lightman (1979), the free-free 
absorption coefficient is 

a ff = 3.7 x 10 8 T- 1/2 Zh^n e z/" 3 (l - e - hu/kT ) g ff , (64) 

where 

^=£^ = ( 65 ) 



op 



In the case of a hydrogen-helium plasma this factor is conveniently represented just by gas 
density, i. e. 

P 



Z 2 rn =n H + ^n H e = — (Y H + Y He ) = — , 



m p m p 



which yields the equation (27). In general, one should use expressions (64) and (65) to find 
the correct form of aff far relevant to the specific chemical composition. 

To be more instructive we conduct such a modification for the case when the plasma 
has a substantial carbon fraction. Using (62), (63) and (65), we write for hydrogen- helium- 
carbon gas 

12 4 
^ = 2AY H + 9Y He + 7Y C = 8-5^-17^/3' (66) 

n^fl-^^V (67) 



m p 



and 



Z 2 ni = J-(Y H + Y He + 3Y C ) = -^(1 + 2Y C ). (68) 



m p m p 



Correspondingly, in all formulae the factor (2 — Y He ) will be replaced by (2 — Y He — Y c ) and 
(8 — 5Y He ) by (8 — 5Y He — 171^/3) • Additionally, the right-hand side of equation (27) has 
to be multiplied by the factor of (1 + 2Yc). Clearly, this modification will add the fifth free 
parameter Yc to the model. Using the general methodology outlined in this paper one should 
be able to produce solutions for the parameter p and the dilution factor. The problem which 
can arise from the inclusion of heavy elements is the possibility for heavy ions to be only 
partly ionized. The ionization degree can also vary throughout the atmosphere. Because full 
ionization and constancy of the gas's chemical composition are the basic assumptions of the 
adopted approach, we cannot explicitly include the effect of ionization in our model. Instead, 
it can be accounted for in a manner similar to our temperature profile correction. First, the 
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approximate atmospheric profiles can be obtained by assuming full ionization. Then, the 
ionization degree can be calculated by solving the Saha equation, and using this solution as a 
zero-order approximation of the atmospheric temperature and the electron number density 
profiles. Finally, one should proceed by solving the hydrodynamic problem, in which the 
partial ionization of heavy elements is taken into account. 

For reasons mentioned above, it is also a problem to include the proper physics for 
the transport of heavy nuclei to the outer layers. Two major processes can contribute to 
this element flow. Bulk motion mixing should dominate in the convection zone close to the 
bottom of the atmosphere. In the outer layers, a strong radiative push should govern the 
process, because of the large resonance cross-sections of the heavy elements. The general 
problem of heavy ions mixing is quite difficult and requires a rigorous approach, which is 
out of scope of this paper. 

As far as the boundary conditions are concerned, modeling of carbon nuclear flashes will 
require higher bottom temperatures. Temperature of the carbon burning zone is argued to 
be about 10 10 K [see Strohmayer & Brown (2001)], which is close to the upper boundary for 
the bottom temperature T b used in our calculations. No peculiarities of the approach were 
detected in the case of very high bottom temperatures. Extremely high temperatures will 
require the correct form of the opacity coefficient k (Paczynski 1983), instead of equation 
(2), which represents a simplified formula for k in the case of modest temperatures. 

Another important issue is the correct accounting for the line emission of heavy elements 
which is detected in the spectral analysis of super-bursts. Strohmayer & Brown (2001) 
argued that this phenomena is due to reflection from the accretion disk during the burst. 
One can estimate the disk heating time by using the standard Shakura-Sunyaev accretion 
disk model (Shakura & Sunyaev 1973) and the fact that approximately 10% of the burst 
luminosity is absorbed by the inner part of the disk (Lapidus & Sunyaev 1985). Simple 
estimates give a timescale of less than a second assuming a mass-accretion rate of the order 
of Eddington or less for the disk accretion regime, and a burst luminosity greater than 5% 
of Eddington, which is detected during several thousand seconds of observation of the super- 
burst in 4U 1820-30 Consequently, the observed spectral feature of the Ka line should rather 
be generated in the burst atmosphere than in the disk. The disk gains the temperature of 
the X-ray radiation very quickly. 

Unfortunately, the origin and behavior of the spectral line features still remain unex- 
plained. The authors plan to inlcude the spectral line effect in the relaxation method, in 
order to calculate the line emission during the X-ray burst and to compare this with the 
observed spectra. 
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Relativistic effects are usually negligible during the strong X-ray burst due to the signif- 
icant radial expansion and the fact that the outgoing spectrum formation occurs at the outer 
layers of the atmosphere. General relativistic effects become important at the contraction 
stage when the extended envelope recedes close to the NS surface [see Lewin, Paradijs & 
Taam (1993)]. Haberl & Titarchuk (1995) applied the full general relativity approach for 
a derivation of NS mass-radius relation in 4U 1820-30 using EXOSAT observations and the 
T94 model. 



7. Conclusion 

This paper follows a common idea of the last decade to fit observational and numerical 
spectra with some model, mostly black-body shapes, to obtain spectral softening/hardening 
factors (London et al. 1986). We improve this technique in several ways. We use for fitting 
more realistic non-blackbody spectral profile, which accounts for the observed power law soft 
excess of X-ray burster spectra. The temperature profile is corrected by solving the tempera- 
ture equation. The existence of the isothermal photosphere during X-ray bursts is confirmed 
numerically and analytically. Finally, we analytically obtain the multiplicative (dilution) 
factor which is not a parameter of fitting anymore but self-consistently incorporated in the 
model. 

We show how the theoretical study of radiatively driven wind phenomenon can produce 
useful techniques for analyzing observational data. It can fulfill the needs of new emerging 
branches of observational X-ray astronomy such as a very promising discovery of super-bursts 
(Strohmayer & Brown 2001), which exhibit photospheric expansion and spectral modifica- 
tions relevant to extended atmospheres. We present the analytical theory of strong X-ray 
bursts, which include effects of Comptonization and free-free absorption. Partly presented 
in some earlier publications, this area of the study of the X-ray burst spectral formation was 
lacking a detailed and self-consistent account. We use numerical simulation to validate our 
analytical theory and to link our solution to energy axes. We show how this information 
can be extracted from spectral data. We provide the analytical expression for the X-ray 
burst spectral shape, which depends only upon input physical parameters of the problem: 
neutron star mass, radius, surface temperature and elemental abundance. Expressions for 
color ratios and dilution coefficient are also given. 

Authors thank Peter Becker for valuable comments and suggestions which improved the 
paper. We are greatful to Menas Kafatos for encouragement and to Center for Earth Science 
and Space Research (GMU) for the support of this research. We appreciate the thorough 
analysis of the presented work by the referee. 
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A. Analytic solution for the radiative transfer problem 



We look for the solution of equation (35) in the form 



J(t,x) = ^— ) B v (T th ) + J(t,x). 



(Al) 



The basic idea is to separate the high-frequency (diluted black body) and the low-frequency 
J(r, x) parts of spectrum, where different physical processes dominate. Kompaneets operator 
L v acting upon black body shape vanishes and we neglect L V {J), which allows us to get the 
solution of radiative transfer problem analytically. At this point r t h is a parameter of the 
problem. The algorithm of determination of r th will be described separately. Substituting 
(Al) into (35) we found for J(t,x) 



d I Id J 
dr I t dr 

with a boundary condition 



3 a // 
r «y 



J = - 



t (Xt 



B„ 



1 - 



4 

TthJ 



T=T t h 



The solution satisfying this condition is presented by 



J(t,x) 



\ r T th 



P w yi{T) J y^ T V^ dT > 



where p(r) = ^ and W{r) is the Wronskian 

W - 

Thus the product 



yi V2 



7 



pW = — 



Functions yi(x) and y-i [x) are 



yi(r)=r/ 4/7 (Jv^^j, 
y 2 (r) = tK 4/7 (^^M(x)t 7 I^ , 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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where I v {x) and K v (x) are modified Bessel function of the first and the second types respec- 
tively. 



The function /(r) in (A4) is the right hand side of equation (A2), namely 
f(r) = --^- = -3r^m!(x)B v 



T (Xt 



8/7' 



eft. 



We introduce a new variable t 
and rewrite solution (A4) as 

rtth 

j(t,x)=B„tV 7 I 4/7 (t) / t 3 / 7 K 4/7 (t) 

Jo 

Using the properties of modified Bessel functions 

J x p K p ^dx = -x p K p + C and K p = K_ p 
we evaluate the integrals in (A10) 

/t t h _ t t h r i-\ 

t^K 4/7 (t) dt = -t^K 3/7 (t) Q = -±LL - tfK 3/7 (t th ), 

£ h t^K A/7 {t) dt = -t^K 11/7 (t)\^= T (^j 2^ - C 7 K u/7 (t th ). 
Finally J(t, x) takes the form 



J(t,x) = B u t 4/7 h/ 7 (t) 



r(f) 3/7 



2 4/7 



8/7 



+t;r(K 11/7 (fe)-%(t tt )) 



(A8) 



(A9) 



(A10) 



J(t,x) = B„t 4/7 h/7(t) 



r(f) r 8 ■ 

"^47? TsTr + 7 A 4/7l^J 



(All) 



The last formula is a solution of equation (35). We can simplify this form by noting that we 
are interested in the solution in the outer layers of atmosphere (emergent spectrum) where 
t — > and t — > 0, and we can use the asymptotic form for small arguments 

Ip{x) w f>Ti) (I) • 
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Making this substitution and putting the result into expression for J(r, x) we find that second 
term in J(t,x) cancels with diluted blackbody term in J(t,x), which takes the form 



J(t,x)=B, 



t 8/7 



r (f) , 8 , -4/7 K ( v 



(A12) 



B. Solution of the temperature equation 

Substituting relation (41) into the equation of radiative diffusion, multiplying it by x 2 
and integrating over the energy range from x* to oo we get 

1 (d 2 R ldR\ f°° x 2 , r „, , . f°° x 2 a ff , 

- -— — / dx — \R(r) - 1 / 11 dx, Bl 

where integrals can be approximated as 

roc 2 roo 

/ ps / x 2 e~ x dx = 2 

J x . e*-l Jo 

and, noting that aff /ax = D^(x) r 3 / 2 ps D r 3 ^ 2 g(x) jx 2 we obtain 

r°° ™2 „ roc 

/ ^^Idx^Dr^ 
Jx, e x -1 a T A, 

Here we used the fact that g(x) ps | ln(2.35/x). The equation for R(r) gets the form 



■^Idx^Dr^ 2 r^dx^ l -\n 22 -^Dr" 2 . 

Ob 4 X J; 



^ - = § In 2 ^D^ 2 [R(r) - 1] = D^ 2 [R{r) - 1]. (B2) 

Boundary conditions for this equation are 

r ^ R(t)^0, 

r^oc i2(r)-l. (B3) 



A general solution of equation (29) is 



R(r) = 1 + rZ 4/7 (^iy/Dr 7 /^ 



(B4) 



where Z i / 7 (z) = ciK 4: /- 7 (z) + 02/4/7(2). In derivation of this formula we take into account 
a well known theorem from ODE theory that general solution of inhomogeneous ODE is 
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the sum of a general solution of the corresponding homogeneous ODE and some particular 
solution of the inhomogeneous ODE, which is chosen equals to unity in our case. The second 
boundary condition and the fact that 

K 4 / 7 (z) ^ z — > oo, 
14/7(2) ^00 z — > 00 

leave only c\ nonzero, and the first boundary condition gives us the value for ci, namely 

1 2 3 / 7 



Cl 



lim tK 4/7 (fv^Dr 7 / 4 ) r (I) W 



-4/7 



where we put r t h = (jj\/~Dj . Then R(r) reduces to 



2 3 / 7 r 

1 7T\ Ka 



t_\ 7/4 ' 

TthJ 



r(|)r«„- 1/7 

Taylor series expansion of K4/7 over r/r^ yields for -R(r) useful relation 



(B5) 



R(t) 



r 



\ 2 



(B6) 



C. Condition at the sonic point (derivation of Y s ) 

We can rewrite the partial derivative in (18) using the obvious relation 

/aP\ dP_ fdT\ dP_ fdy\ 
\dp) s ~ df \dp) E + ~d^ \dp) s 

Differentiating the equation of state (8) we obtain derivatives of pressure 

OP / 1\ 4aT 3 , OP 1 aT 4 
— — = I 1 — I — — J and — = 



dT \ yj 3 dy y 2 3 ' 

and differentiating (21) with respect to p we get 

(dT\ = _T(l_ + A \(dy 



\dp J 3 A \y / \dp J E 
From the other hand differentiation of (9) gives us 

apm p f3T 2 dT T 3 dy 



3k \ y dp y 2 dp 



I. 
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Combination with the previous equation it yields: 



dy\ 3k y 1 

dp) E ~ ~ap,m p T 3 \\ + 3(1 + Ay) J ' 



3A; 



2/(1 + 4y) 
\dpj s apm p T 2 \X + 3(1 + Ay) 



Now, combining all found derivatives we have 



k 



(dP_ 

\dp J 3 \im v 
or in a more compact form 



4 1 + 



y{l + Ay) 



yj \\ + 3(1 + Ay) J A + 3(1 + Ay) 



A 



T 



k 



dp J s jj,m p 



A + 4(1 + y)(l + 4y) 



A + 3(l + 4y) 



T, 



which is, in fact, a sonic point condition (18). 



(CI) 



D. Reduction of Hydrodynamical Problem to First-Order ODE 

We derive an expression for the derivative of velocity v y through v and y. 

We substitute the temperature profile found in (21) to (5) to obtain p as a function of y 



p = p{y) 



y-s/*- 1 exp 

3A; 



12(2/ -1) 
A 



Using this expression for p(y) and equation (12) we get 



1/2 



(■ 



3§k 



Then we get derivatives 



1/2 



t/2A^2 exp 



6(y -i; 

A 



_ i 
f 2 . 



3 1\ 1 6 
2A + 2 J y + A 



(if 2t> 
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Differentiation of (21) also gives us 



dy 



cxp 



A 



1 



By a combination of all these derivatives we obtain 



dT _dTdy _dT (dr dr dv 
dr dy dr dy \dy dv dy 

Substitution of it into (5) yields 

lQnckr 2 y dT 



-i 



2T(% + 1) 
r\y 



3 \ 1 12 v r 
A ) y X 



32nck 



(Ay + 1)(1 + aT)Tr 



lim p K dr iim p \K (2 — Y He ) [(3/A + 1) 1/y + 12/A — v'/v] 
But from (13) we also have 

L r = *-$ [h+- " 

\ 2 r 

Equating the last two expressions for L r we finally find v y as 

v' -f(vv)-v\(l + 3 I±^ 1 - 75 2 ^(8-5y ge )(l + 4 y )(l + aT) 

Here and T^g represent neutron star radius and the bottom temperature of atmosphere 
in terms of 10 6 cm and 10 9 K correspondingly. 
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Table 1. Parameters for model m = 1.5, R ns = 13.5 km, Y He = 1 





kT 


T/T 


Tth 


P 


<j> a 




T s 


v s 


r s 


r p h 


10 9 K 


keV 










10 22 g 


0.1 keV 


10 3 km/s 


10 3 km 


10 3 km 


7.0 


0.197 


0.099 


45.1 


1.56 


93.9 


173.3 


0.21 


1.31 


57.9 


4.16 


6.5 


0.208 


0.105 


44.7 


1.64 


87.2 


128.9 


0.22 


1.39 


51.6 


3.71 


6.0 


0.221 


0.111 


44.1 


1.75 


80.5 


93.6 


0.24 


1.48 


45.6 


3.28 


5.5 


0.236 


0.118 


43.3 


1.87 


73.8 


66.1 


0.25 


1.58 


39.9 


2.88 


5.0 


0.252 


0.127 


42.4 


2.02 


67.1 


45.1 


0.27 


1.70 


34.5 


2.50 


4.5 


0.272 


0.137 


41.5 


2.20 


60.4 


29.6 


0.29 


1.84 


29.4 


2.14 


4.0 


0.295 


0.148 


40.3 


2.42 


53.7 


18.5 


0.32 


2.01 


24.6 


1.80 


3.5 


0.324 


0.163 


39.0 


2.71 


46.9 


10.9 


0.35 


2.22 


20.1 


1.48 


3.0 


0.359 


0.180 


37.4 


3.09 


40.2 


5.86 


0.39 


2.50 


16.0 


1.19 


2.5 


0.406 


0.204 


35.6 


3.60 


33.5 


2.83 


0.44 


2.85 


12.2 


0.92 


2.0 


0.470 


0.236 


33.4 


4.34 


26.8 


1.16 


0.51 


3.36 


8.83 


0.67 


1.75 


0.510 


0.257 


32.2 


4.85 


23.5 


0.68 


0.55 


3.69 


7.29 


0.56 


1.5 


0.565 


0.284 


30.8 


5.53 


20.1 


0.37 


0.61 


4.12 


5.86 


0.45 


1.25 


0.634 


0.318 


29.2 


6.44 


16.8 


0.179 


0.68 


4.69 


4.53 


0.36 


1.1 


0.686 


0.344 


28.0 


7.18 


14.8 


0.108 


0.74 


5.12 


3.80 


0.30 


1.0 


0.727 


0.365 


27.2 


7.78 


13.4 


0.074 


0.78 


5.47 


3.33 


0.27 


0.9 


0.775 


0.389 


26.3 


8.51 


12.1 


0.049 


0.84 


5.87 


2.89 


0.23 


0.8 


0.831 


0.417 


25.3 


9.40 


10.7 


0.031 


0.90 


6.36 


2.46 


0.20 


0.7 


0.898 


0.451 


24.3 


10.5 


9.4 


0.018 


0.97 


6.95 


2.06 


0.17 


0.6 


0.981 


0.492 


23.0 


12.0 


8.0 


0.010 


1.06 


7.69 


1.68 


0.14 


0.5 


1.086 


0.545 


21.6 


14.0 


6.7 


0.005 


1.17 


8.65 


1.33 


0.11 


0.4 


1.224 


0.614 


19.9 


17.0 


5.4 


0.002 


1.32 


9.95 


1.01 


0.09 


0.3 


1.417 


0.711 


17.8 


21.9 


4.0 


0.001 


1.53 


11.8 


0.71 


0.07 



a <I> is in units of the critical mass-loss rate, i. c. divided by L_b/c 2 
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Table 2. Parameters for model m = 2.6, it! ns = 14.0 km, Y He = 1. 





kT 


T/T 


Tth 


P 


$ 




T s 


Vs 


r s 


r p h 


10 9 K 


keV 










10 22 g 


0.1 keV 


10 3 km/s 


10 3 km 


10 3 km 


7.0 


0.248 


0.110 


45.1 


2.12 


56.2 


115.8 


0.242 


1.80 


53.4 


3.53 


6.5 


0.261 


0.116 


44.5 


2.24 


52.2 


86.1 


0.256 


1.90 


47.7 


3.15 


6.0 


0.277 


0.123 


43.6 


2.40 


48.2 


62.6 


0.271 


2.02 


42.2 


2.80 


5.5 


0.294 


0.131 


42.6 


2.58 


44.1 


44.2 


0.288 


2.16 


36.9 


2.46 


5.0 


0.313 


0.139 


41.5 


2.80 


40.1 


30.2 


0.308 


2.32 


31.9 


2.15 


4.5 


0.337 


0.150 


40.3 


3.07 


36.1 


19.8 


0.331 


2.52 


27.2 


1.84 


4.0 


0.364 


0.162 


39.0 


3.40 


32.1 


12.4 


0.359 


2.75 


22.8 


1.56 


3.5 


0.398 


0.177 


37.4 


3.82 


28.1 


7.27 


0.394 


3.04 


18.7 


1.29 


3.0 


0.440 


0.196 


35.7 


4.36 


24.1 


3.93 


0.437 


3.41 


14.9 


1.04 


2.5 


0.495 


0.220 


33.8 


5.12 


20.1 


1.90 


0.494 


3.89 


11.4 


0.81 


2.0 


0.570 


0.254 


31.5 


6.21 


16.1 


0.78 


0.572 


4.58 


8.24 


0.59 


1.75 


0.620 


0.276 


30.2 


6.97 


14.0 


0.46 


0.623 


5.03 


6.80 


0.50 


1.5 


0.682 


0.304 


28.7 


7.97 


12.0 


0.25 


0.687 


5.61 


5.47 


0.40 


1.25 


0.763 


0.339 


27.0 


9.33 


10.0 


0.121 


0.770 


6.38 


4.24 


0.32 


1.1 


0.824 


0.366 


25.9 


10.4 


8.83 


0.073 


0.833 


6.96 


3.56 


0.27 


1.0 


0.872 


0.388 


25.1 


11.3 


8.03 


0.050 


0.883 


7.43 


3.12 


0.24 


0.9 


0.928 


0.413 


24.2 


12.4 


7.22 


0.033 


0.940 


7.98 


2.71 


0.21 


0.8 


0.993 


0.442 


23.3 


13.7 


6.42 


0.021 


1.008 


8.63 


2.31 


0.18 


0.7 


1.072 


0.477 


22.2 


15.4 


5.62 


0.012 


1.089 


9.43 


1.94 


0.15 


0.6 


1.169 


0.520 


21.0 


17.6 


4.82 


0.007 


1.188 


10.4 


1.59 


0.13 


0.5 


1.292 


0.575 


19.7 


20.6 


4.01 


0.003 


1.313 


11.7 


1.26 


0.10 


0.4 


1.455 


0.647 


18.1 


24.9 


3.21 


0.001 


1.479 


13.5 


0.95 


0.08 


0.3 


1.682 


0.749 


16.2 


32.0 


2.41 


0.0005 


1.710 


16.0 


0.67 


0.06 



